Multi-phase matching in the Grover algorithm 
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Phase matching has been studied for the Grover algorithm as a way of enhancing the efficiency of 
the quantum search. Recently Li and Li found that a particular form of phase matching yields, with 
a single Grover operation, a success probability greater than 25/27 for finding the equal-amplitude 
superposition of marked states when the fraction of the marked states stored in a database state is 
greater than 1/3. Although this single operation eliminates the oscillations of the success probability 
that occur with multiple Grover operations, the latter oscillations reappear with multiple iterations 
of Li and Li's phase matching. In this paper we introduce a multi-phase matching subject to a 
certain matching rule by which we can obtain a multiple Grover operation that with only a few 
iterations yields a success probability that is almost constant and unity over a wide range of the 
fraction of marked items. As an example we show that a multi-phase operation with six iterations 
yields a success probability between 99.8% and 100% for a fraction of marked states of 1/10 or 
larger. 

PACS numbers: 02.70.-c, 03.65.-w, 03.67.Ac, 03.67.Lx 
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I. INTRODUCTION 



The quantum search algorithm introduced by 
Grover [J 0, d, 01 constitutes a major advance in 
quantum computing. It enables us to find a marked 
state stored in a database state consisting of N un- 
ordered basis states in only 0{^/N) Grover operations. 
A number of modifications and generalizations of 
the original Grover search algorithm have been pro- 
posed p, Hi 0, B B [13, EH ■ In particular, phase matching 
methods in the Grover algorithm have been extensively 
examined 0, 0, The outcome of the search algo- 

rithm is characterized in terms of P{X), the probability 
of obtaining an equal-amplitude superposition of the 
marked states where A is the ratio of the marked states 
to all the states stored in the original database state. 

Recently, Li and Li (lo| proposed a new phase match- 
ing for the Grover algorithm and they obtained an im- 
proved success probability P(A) over a wide range of the 
ratio A. They introduced the set of the Grover operators 
(details are described in Eqs. (P?T|) and ((^ ): U = I - 
(l_e")^f^-i \ti){ti\ 'AndV = /e''^ + (l-e*'5)|0®")(0®"|. 
The phase factor e*^ in the first term of the operator V 
was first introduced in Ref. [To| . In the new phase match- 
ing the number of phases is the same as the usual one but 
the form of the phase shift operator V is different. Li and 
Li found the remarkable result that a single Grover oper- 
ation of the new phase matching yields P{X) > 25/27 for 
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1/3 < A < 1. This is significant in the sense that with 
only one Grover operation the efficiency of the Grover 
algorithm is substantially improved in the range of val- 
ues of A where the efficiency of the original algorithm 
deteriorates. 

This phase matching has another interestin g a spect 
that was not explicitly pointed out by Li and Li [l3| . For 
a given values of A in the range 1/4 < A < 1 , one Grover 
operation with the phases a = —(3 = arccos(l — 1/2A) 
yields exactly P = 1. [See Eq. (|2.11[) in the following.] 
This results was obtained earlier by Chi and Kim [l^l who 
considered a modified Grover operator of arbitrary phase. 
The special case of A = 1/2 yields a = — /? = ±7r/2, 
which are the phases found in Ref. [l^l- This aspect of 
the phase matching is also significant because it implies 
that one can always find the equal-amplitude superposi- 
tion of the marked states by only one Grover operation 
when A is greater than 1/4 by tuning the phases a and 
(3 appropriately for the given A. Conditions for a suc- 
cess probability of unity have been studied by previous 
authors. See, for example, Refs. [Tsl. [T3|. 

It should be pointed out, however, that the so-called 
new phase matching of Ref. [l^l is equivalent to the orig- 
inal phase matching of Long et al. [y] . When the second 
operator is defined as V = e~'^^V, it becomes the phase- 
matching operator of Long et al. The only difference 
between the two is that the overall state is multiplied by 
a phase factor and so the amplitudes of the components 
are different, but the probabilities are the same. Thus 
the remarkable result of Li and Li can also be seen to 
follow from the operator of Long et al. Analytically the 
formulation by Li and Li is somewhat more transparent 
and hence we use it throughout this paper, except in the 
Appendix where we explicitly show the equivalence of the 



two formulations by calculating the probability profile. 

Thus a number of aspects of the Grover algorithm with 
phase matching, already alluded to, are of particular in- 
terest and they form the objectives of this study. We fo- 
cus on high success probabilities with as few iterations as 
possible in order to enhance the efficiency of the quantum 
search. We emphasize the following three objectives: (1) 
the elucidation of features of the phase-matched Grover 
operations with a small number of iterations that yield 
success probabilities P(A) close to one over a wide range 
of values of A, (2) given a value of A the determination 
of the phase-matched Grover operator(s) that results in 
P(A) = 1 exactly, and (3) the elucidation of the features 
of the phase-matched Grover operators that allow us to 
obtain P(A) = 1 for very small values of A. 

In this paper we explore the search algorithm with 
these objectives in mind using the advantages of a few 
multiple Grover operations with phase matching. It is 
well known that a multiple application of the original 
Grover operation gives rise to intensive oscillations of P 
as a function of A and such oscillations deteriorate the 
efficiency of the algorithm. This undesirable feature re- 
mains even in the new phase matching of Li and Li, as 
we will illustrate. We show that if we introduce a multi- 
phase matching subject to a certain matching rule, we 
can obtain a multiple Grover operation that yields a suc- 
cess probability almost constant and unity over a wide 
range of A, e.g., 0.1 < A < 1. This is also significant in 
the sense that when A is greater than a small minimum 
value we can always find the superposition of the marked 
states with high degree of certainty without (re)tuning 
the phases. 

In the next section we set up the algorithm of the multi- 
phase matching in the framework of the phase matching 
of Li and Li [lO| and analyze the efficiency of the algo- 
rithm by considering a single matched phase and a two- 
stage multi-phase matching. We also obtain an exemplar 
of a good probability profile for a six-stage multi-phase 
matched operator. In Sec. IIIII we consider the success 
probability for small A by using the Grover operations 
with a phase other than it. We summarize our results in 
SecHVl 

II. MULTI-PHASE MATCHING IN THE 
FRAMEWORK OF THE NEW PHASE 
MATCHING 

The new phase matching in the Grover algorithm pro- 
posed by Li and Li [lo[ is defined with the two operators, 

M-l 

f/ = /-(l-e-)^K,)(i,| (2.1) 
;=o 

I 



y = /e''' + (l-e^'3)|0®")(0®"|. (2.2) 



where jO®") is the n-qubits initial state, M is the num- 
ber of target (marked) states stored in an unstructured 
database state, and the \ti) denote the target or marked 
states. The database state is given as {cj)) = ij'»"|0'*"), 
where H is the Walsh-Hadamard transformation. The 
state \4>) is an equally- weighted superposition of the 
TV = 2" basis states, \wi), l = 0,...,N-l. The fraction 
A of the target states is defined as A = M/N . The U 
and V of Eqs. (12.11) and (|2.2p are both unitary as was 
shown in Ref. [10|. With a = /3 = tt, J7 and V re- 
duce to the Grover operators of the original algorithm. 
As we mentioned in Sec. U Li and Li showed explicitly 
that a single Grover operation of the new phase match- 
ing (ii-»»yiy®»){777®"|0®") with a = -/3 = tt/2 yields 
a success probability P(A) > 25/27 for 1/3 < A < 1. 

We introduce a multi-phase matching within the 
framework of the new phase matching. We rewrite the 
database state |0) = i/®»|0«5») = A^"^/^ E/loM^i) in 
terms of A as 




= \/T^|i?) + \/A|T), (2.3) 



where 



^ N-M-l M-l 

\R) = , y In), \T) = ^y \ti). 

(2.4) 

The state \T) is the uniform superposition of the marked 
states and \B) is that of the remaining states |r/). 
They are both normalized to unity and orthogonal 
to each other. In the following, for convenience, we 
work in the two-dimensional space defined by the ba- 
sis |r)}. The two-dimensional representations of U 
and H®''VH®'' = le't^ + (1 - e*'')|(/))(0| are 



U 



1 

e" 



(l-e'^)(l -A)-fe'^ (l-e'/5)v/A(l-A) 
{l~e"^)y^X{l- \) {l~e'^)X + e"^ 



(2.5) 
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We write the multiple Grover operation with the multiple phases aj and Pj (j = 1, . . . , fc) as 



2j =G(afc,/3fc)GK-i,/?fc-i)---G(ai,/3i) (\^) , (2-6) 
where one Grover operation G(aj, (3j) (j = I, k) in this representation is 



The success probability of finding the superposition of 
target states is given by Pfe(A) = \dk\'^. 

We now consider the one- and two-pair-phase cases 
before increasing the phase-matching to six different pairs 
of phases in order to obtain P(A) nearly equal to unity 
over a large range of values of A. In other words, we 
discuss the k = 1 and the k = 2 cases in detail first, and 
then proceed to the numerical results of the fc = 6 case. 



A. Multi-phase matching with one pair of phases 



When fc = 1 , Eq. (j2.6|l reduces to 



(2.8) 



Since A lies between zero and one, the range of a is 7r/3 < 
a < TT. The boundary point of this range a = 7r/3 occurs 
when P(A = 1) = 1, and similarly a = tt when P(A = 
1/4) = 1. 

We can express P(A) as a function of A depending on 
the parameter a, 

P(A) = 

= A [5 — 4 cos a — 4 (1 — cos a) (2 — cos a) A 

+4(l-cosa)2A2]. (2.12) 

For a = 7r/2 the equation reduces to Eq. (14) of Li and 
Li [l^l . Since Eq. (|2.12p is cubic in A we expect a local 
maximum and a local minimum in the range < A < 1 
at Aniax and Amin respectively, where 



Since we first focus on cases of complete success P = 
= 1 — = 1. we can equivalently consider the 
condition ui = 0. In general 



A. 



1 



5 — 4 cos a 



Ul 



Vr^[l-A + e'^^^A+(e--e'^(-«)A]. (2.9) p(a_) . i, P(A_) - ^ 



2(1 — cos a) 6(1 — cos a) 

Furthermore the extrema are 

1 + cos a) (5 — 4cosq;)^ 



(2.13) 



The condition that ui be zero leads to 

— = 1 — cos a — cos f3 — cos(q! + /3) 
A 

+i[sin (a + /3) — sina — sin/3]. (2.10) 

The fact that A must be real implies that (1) f3 = —a, 
(2) either a or /3 are zero, or (3) both a and (3 are zero. 
When /3 = in the operator V of Eq. (|2.2p , the operator 
is the identity and the overall effect of operator U of 
Eq. (|2.ip by itself would cause the phase of the marked 
states to be changed, but the probabilities of marked and 
unmarked states would remain the same. When a = 0, 
then U = I and G = H®''VH®'\ The initial state \(j)) 
is an eigenvector of G with eigenvalue 1. Thus G does 
not cause any evolution in |^). The success probability 
is P = A, which is the success probability of the classical 
algorithm. As no quantum improvement to the search 
algorithm is achieved, we eliminate the case of a = and 
any nonzero /3 from the solutions of Eq. (|2.10p . Thus 
only the solution /3 = — a is meaningful, and yields, as 
mentioned in Sec. [T] and in Ref. [12], P = 1 when 



27 (1 — cos a) 



(2.14) 

We illustrate different cases in Fig. [TJ It is evident that 



P(k) 




a=7t/3 
a=it/2 
a=3i/2+0.1 
2a=it/3 
5a=7t/6 
a=jt 



0.4 0.6 



0.8 



FIG. 1: (Color online) Plots of P(A) for different values of the 
parameter a for the case of one iteration. (See Eq. (I2.12|l .) 



-/? = arccos(l - 1/2A). 



(2.11) the a — it /2 case, which is the one used by Li and Li [T^], 
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gives the optimal profile for the success probability. Op- 
timal here could be defined as the largest average P over 
the range of A, or the largest range of A over which 
P > 25/27. 

B. Multi-phase matching with two pairs of phases 

We now consider Eq. (|2.6p for k — 2 and we again con- 
centrate on the upper component of the vector (w2, ^2)"^- 
The general expression for it is too lengthy to give here, 
but again we demand that for an arbitrary value of A the 
imaginary part is zero to obtain the matching relation- 
ship for the phases. Apart from a factor of \/l — A the 
expression of Im U2 contains a term in A and another in 
A^. Demanding that the coefficients of each power of A 
vanishes gives us two equations involving ai, 0.2, Pi, and 
(32- Solving for (3i and P2 in terms of ai and 02 we obtain 
the following four solutions: 

{Pi = -ai, P2 = 0} 

= -a2, P2 = -ai} 
{Pi = ^2 - 0} 
{Pi = 0, /32 = -(ai+a2)}. 

Since one of Pi and P2 is zero for the first and last so- 
lution, the operation is then reduced to one iteration. 



and for the third solution the two iterations would not 
change the probabilities of the marked and unmarked 
states. Thus the only solution that gives new informa- 
tion is the one where Pi = —a2 and P2 = —ai. (The 
fact that ImM2 = is a necessary, but not a sufficient, 
condition for this solution.) After obtaining the matched 
phases for which ImM2 = 0, we set Rc'U2 = to solve for 
the values of A which gives P = 1. 

The expression for U2 is then real and can be written 

as 

U2 = {1 + 2 [(1 — cosQ!i)(— 2 + cosq;2) ^ sinai sina2] A 
+4 (1 — cosai)(l — cosa2 )A2}V1~A. (2.15) 



The factor multiplying vT— A is quadratic in A and 
hence it can vanish for two values of A. Thus we can 
ask ourselves the questions, suppose two values of A be- 
tween zero and one are given at which P{X) = 1, what 
are the corresponding values of a and what limits are 
there on the possible values of A that satisfy P(A) = 1? 
If Ai and A2 are the roots of the equation 

u2iX)/VT^=0, (2.16) 

then cos ai and cos a2 satisfy the equations 



SAiAaCOS^as + [4(Ai + A2)(l - Ai - A2) - 8AiA2]cos2 + [8(Ai + A2)^ - 12(Ai -I- A2) - 8A1A2 +4] cos as 



-4(Ai -t- \2f + 8(Ai -h A2) - 5 + 8A1A2 = 0, 



(2.17) 



COSCki = 1 



1 



4(1 — cosq;2)AiA2 ' 



(2.18) 



In order to have a sense of the values of ai and a2 that are 
valid, we have minimally the condition that the discrimi- 
nant of Eq. (j2.16p (quadratic in A) should be nonnegativc 
to avoid complex values of A. In Fig. [2] we plot the dis- 
criminant as a surface z = D{ai, 012); the intersection of 
the surface with the xy plane gives the boundary of the 
non-allowed ai and 02 values. 



we obtain a change in the smaller A of 



Given Ai and A2 one can solve Eq. (|2.17p for cosq!2 
and using it we obtain cos ai from the second equation. 
Only those solutions that yield real angles ai and a2 
arc meaningful for the unitary operators. The minimum 
value of A for which P = 1 occurs when ai = a2 = tt. In 
that case A = (3-V5)/8 = 0.09549. It can be shown that 
varying ai or a2 by a small amount away from tt always 
leads to an increase in the A which corresponds to the 
smaller of the two values of A. When we let ai 2 = Tr+ei 2 



AA 



1 
160 



2V5ei 



5 - 3\/5 



62 +{22-8V5)el 



V2VE 

(2.19) 

which is positive regardless of the signs of £1,2- The larger 
A can increase or decrease with changes in the phases(s). 

We obtain a particular example using the procedure 
described above. We search through combinations of A1.2 
and find that Ai = 2/5 and A2 = 4/5 give good results. 
In this case ai = 1.00889485 and 02 = 2.30794928. We 
find local minima of P(A) at A = 0.5767 and A = 0.9433 
at which P = 0.9936 and 0.9966, respectively. The cor- 
responding graph of the success probability as a function 
of A obtained with the two-stage multi-phase operator is 
shown in Fig. [3] and compared with double iterations of 
the Grovcr operation and that of Li and Li [l3| . 

It would be interesting to examine a classical counter- 
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FIG. 2: (Color online) The discriminant of Eq. (|2.16p as a 
function of Ofi , and 02 ■ The curve of intersection of the surface 
with the D = plane separates the smaller values of Qi for 
which there are no physical operators. 



part of P(A). The probability of failing to find one of M 
marked objects out of N objects is {N — M)/N = 1 — A. 
The probability of failing twice in a row is 



(1-A) 



N -l-M 



(1-A) 1 



A 



iv-1 J "' y i-i/N 

The probability of failing k times in a row is 
A \ A 



(1-A) 1 



n 



1-A 1 



1 - 1/A^ 

n - 1 



1 



1 - (/c - 1)IN 



N 



Thus the probability of finding at least one of the M 
items in k successive attempts is 



-Pclassical(A) — 1 — 



n=l 



1-A 1- 



n - 1 
N 



(2.20) 



If ^ TV, this probability is approximately Pciassicai(A) ~ 
1 — (1 — A)*"', which we interpret as the classical counter- 
part of P{\)- This probability with fc = 2 is also plotted 
in Fig. [31 



C. Multi-phase matching with six pairs of phases 

We show that if we match the multi-phase aj and /3j 
(j = 1, . . . , fc) with fc = 6 in accordance with a certain 
matching rule (best fit), we can obtain a multiple Grover 
operation that yields P{\) « 1 in a wide range of A. We 
found this best solution for six Grover iterations by a 
nonlinear fitting to the ideal probability curve P{\) = 1 
for < A < 1- The phases aj and f3j found in this way 
are given in the left side of Table HI 
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FIG. 3: (Color online) Plots of P(A) obtained after two it- 
erations of the multiphase case with qi = 1.00889485 and 
a2 = 2.30794928. For comparison the double Grover iter- 
ation (a.\ — Qf2 = 71"), the double iteration of Li and Li 
{ai = a2 — 7i"/2), and the classical counterpart are shown 
as well. 



j 


aj/iv 


ft 




r>/\ (local min) \ 

-^l-^., ) 


1 


1.20560132 




0.10777 


0.9980 


2 


1.29806396 


-as 


0.23793 


0.9993 


3 


1.31701508 


-04 


0.41889 


0.9996 


4 


1.33356767 


-a-s 


0.62393 


0.9997 


5 


0.47289426 


—a2 


0.81366 


0.9997 


6 


1.66668634 




0.94483 


0.9995 



TABLE I: Phase parameters for the six-parameter multi- 
phase matching, and results for local maxima (P(A) — 1) 
and local minima of P{X)- 



It is remarkable that aj and f3j are matched to each 
other such that aj = — /?6-j+i . The signs of aj and (3j are 
opposite to each other, which is consistent with the case 
of the new phase matching of Ref. [l3|, i.e., the k ~ 1 case 
with ai = —Pi = tt/2. The matching rule aj = —(3k-j+i 
between the multi-phases aj and (3j holds for any k in 
the best solution obtained by the nonlinear fitting to the 
ideal probability curve P{X) = 1 for < A < 1, although 
we omit to show cases other than those for which fc = 1, 
2, and 6. 

Fig. m shows the success probabilities obtained by 
six Grover operations with the multi-phase matching of 
Eq. (|2.6p . The inset of Fig. 2] shows that there are six 
values of Ai at which P{\i) ~ 1 exactly. They are given 
in the right side of Table [J along with local minimum 
values of the function P(A) which occur between 0.1 and 
1. 

We studied the k = 5 case in the same way and ob- 
tained a graph similar to Fig. [4] with P(A) = 1 for five 
values of A other than unity. The local minima of P{X) 
are lower and the minimum value of A for which P(A) = 1 
is slightly larger than in the fc = 6 case. The matching 
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X 

FIG. 4: (Color online) Success probability curves P(A) ob- 
tained with the multi-phase matching with six iteration in 
the operation. The phase parameters are given in the left 
side of Table HI 

rule aj = —f3k-j+i is also satisfied for the fc = 5 case as 
it was for fc = 1, 2, and 6 cases. We are confident that 
for any fc > 1 this matching rule for the best fit holds so 
that in general one finds k values of A for which P{X) = 1 
and the smallest A for which P(A) = 1 decreases as k in- 
creases. 




FIG. 5: (Color online) Success probability curves Pi (A) at the 
ith stage obtained with the six-stage multi-phase matching. 
The phase parameters are given in the left side of Table U 

Returning to the six-stage multiple phase operation, 
we define Pj{X) {j — 1, ... ,6) as the success probability 
curves after j steps of the six-stage multi-phase opera- 
tion. As seen in the Figs. [4] and EJ Pei^) w 1 is achieved 
for 0.1 < A < 1 in the sixth Grover operation. This 
is significant in the sense that if A is greater than 0.1, 
we can always find the superposition of marked states 
by just six Grover operations. In contrast to the shape 
of the curve for Pe{X), Fig. [5] shows that each Pj{X) for 



j = 1, . . . , 5 depends strongly on A and is far from the de- 
sired success probability Pe{X). The curves do not mono- 
tonically approach the desired success probability Pe^X) 
when A > 0.05. In particular, P5{X) is quite different 
from the desired probability PeiX). However, in the final 
(sixth) step the desired probability PeiX) is obtained. 
This is in contrast to the fixed-point iteration schemes 
studied in Refs. Sim, HI. 




X 

FIG. 6: (Color online) Success probability curves -Pi (A), 
P3(A), and P6(A) obtained with the single phase matching 
with aj — —Pj = 7r/2 (j = 1, • ■ ■ , 6). The Pi(A) indicates the 
success probability of the new phase matching of Ref. [lol |. 
The yellow curve is the success probability of the original 
Grover's algorithm, where the graph segments are plotted for 
optimal iteration times indicated by k. Iteration times for 
small A are omitted. 

Figure[6]shows the success probabilities obtained by six 
Grover operations with the single phase matching with 
Uj = —(3j = 7r/2 (j = 1, . . . , 6) , where we showed only 
Pi (A), PsiX) and Pe{X). The Pi (A) is the success prob- 
ability of the new phase matching obtained by Li and 
Li [lal. As stressed in Ref. [3, the success probability 
is substantially improved in A > 1/3 by a single Grover 
operation, compared with that of the original Grover al- 
gorithm indicated by the yellow line, where the probabil- 
ity is plotted for optimal iteration times indicated by fc. 
However, P3(A) and P6(A) obtained by multiple Grover 
operations with the single phase matching show intensive 
oscillations with A. As we have shown, such undesirable 
oscillations can be eliminated by the multi-phase match- 
ing subject to the matching rule aj = — /Jg-j+i. 

Here we should note that the nonlinear fitting is not 
unique. The phases Uj and Pj given in Table U were ob- 
tained by minimizing the function Xi where the Xi ^re 
the differences at A = A^ of the ideal probability P(A) = 1 
and the probability function P6(A). If we take, for exam- 
ple, a function such as \xi\ we obtain another solution. 
Although this solution gives almost the same P6(A), the 
probability curve is shifted slightly toward larger values 
of A, so that the local extrema are also slightly moved to 
the right. Since we emphasize obtaining P(A) w 1 over 



7 



as wide a range of A as possible we adopted the solution 
that uses the Xi for the fitting. 



III. ITERATION OF GROVER'S OPERATION 
WITH PHASE OTHER THAN tt 

In this scetion we consider the repeated application of 
Grover's original operation generalized to have a phase 
other than tt. We focus in particular on cases with small 
A for which the success probability with the multi-phase 
matching is small, and determine the conditions that 
yield success probabilities close to unity. 

Consider a single Grover operation with matched 
phase, Eq. (|2.7p . but with (3 = —a, 



g ^ [ (1 - )(1 - A) + e-' " (e» - 1)VA(1-A)\ 
' 1^ (1-e— )v/A(l-A) (e"-l)A + l )' 

(3.1) 

Note that detGi = 1. We obtain eigenvalues a of the 
matrix Gi by solving 

/(cr) = det(Gi - al) = 0. (3.2) 

The characteristic function /(cr) is 

f{a) =(7^ + 2[-l + (1 - eosa)A]cr + 1. (3.3) 

The equation /(cr) = yields solutions 



cr = 1 — (1 — cos a)\ ± iy^ (1 — cosci!)A[2 — (1 — cosq;)A] 

(3.4) 

We define x as 



X = (1 — cosa)A, 



so that the eigenvalues can be written as 



arctan 



(3.6) 



We choose the definition of the arc tangent so that as x 
varies from to 2, goes from to tt. We can rewrite 
the function /(cr) as 



/(cr) = cr^ - 2crcoS(/) + 1. 



(3.7) 



By the Cayley-Hamilton theorem [T^, page 91] /(Gi) = 
0, so that we obtain the identity 



Gi = 2Gi coscj)-!. 



(3.8) 



This means that Gi iterated any number of times can 
be written as a linear expression of Gi. In fact for k 
iterations it can be shown by induction [l3| that 



G'^ = ^— [Gi sm(k(j)) - sin((/fc - 1)0)1 . (3.9) 
sm (h 



Consider now the k iterations of the Grover operation, 
so that 



Gf 



k I 



Va 



(3.10) 



(3.5) This yields 



Jk. 



sin {k(j)) 



(l-e-*")(l-A) + e-*" (e"-l)v/A(l-A) 
^ (l-e-»)v/A(l-A) (e--l)A + l , 



- sin ((/s- 1)0) / 



r 



(3.11) 



Thus the expression for Uk is 



Uk 



{sin (fc0)(l - 2x) - sin {{k - 1)0)} . 



(3.12) 

We require Uk — so that P = 1. A trivial solution is 
A = 1. We also note that = yields x = — l/(2fc). 
Since x must be positive sin0 7^ 0. Thus we need to 
solve only 



sin (fc0)(l - 2x) - sin ((fc - 1)0) = 0. 



(3.13) 



The solutions are values of a; — (1 — cos a) A for which 
P ~ I. Thus we have P = 1 for combinations of a and 
A. For instance, when a = tt, then A = x/2. In Fig. [3 we 



display the P(A) curves for six (fc — 6) iterations when a 
has different values. 

For large k we can estimate the smallest value of A for 
which P is unity. We rewrite Eq. (|3.13p so that 



tan(fc0) 



cos 6 — 1 + 2x 



(3.14) 



The value of x which is the solution occurs for the x 
coordinate of the point of intersection of the curves rep- 
resented by the left side and the right side of Eq. p.l4p . 
The curve on the right is a smoothly decreasing positive 
function starting at infinity when x = and asymptoti- 
cally approaching the positive a; axis. The curve on the 
left starts at zero and increases to positive infinity when 
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FIG. 7: (Color online) Plots of P(A) obtained after six itera- 
tions with a single a for the cases of a = tt/S, 7r/4, 7r/2, and 

TT. 



=- 11/2. When k is large this occurs for small val- 
ues of 4> or small values of x. Thus using the condition 
(j) Ki 7r/(2fc), we obtain 



■v/ x{2 ~ x) ^ TT 

arctan-^^ ^ — or 

1 - -T 2k 



^^^^<^. (3.15) 

l-x 2k ^ ' 



\~x 

This leads to the approximation of the smallest value of 
X for which P is one as Xinin^7''^/(8/s^); for a = tt (the 
Grover case) Amin = a:niin/2. This approximation leads 
to Aniin = 0.017 for A: = 6 and a = tt, whereas the exact 
solution of Eq. (|3.14[) gives 0.014. As k gets larger the 
approximation improves further. 

A. The Grover algorithm as a special case 

We recover the Grover algorithm starting with 
Eq. (|3.12p and setting a = tt or a; = 2A. Then it 
follows from Eq. (|3.6p that sin0 = ■\/4A(l — A) and 
cos(/) = (1 — 2A). The Ufe of Eq. (j3.12|l can be reduced to 

Wfc = — Va sin(fc(/)) -I- Vl — \ cos(k(^). (3.16) 

Define sin = VA. Then 

life = cos(fc0-t-6'). (3.17) 

We can show that ~ 29, so that 

Uk ~ cos [{2k + 1)9] = cos[(2fc + 1) arcsin(VA)]. (3.18) 

This is Eq. (6) of Rcf. |10|]. Furthermore 

P=l-ul = sin2[(2fc + 1) arcsin(\/A)]. (3.19) 

Thus each iteration effectively rotates the state through 
an angle of 9/2 = arcsin(-\/A)/2. We can use this to 



estimate the number of iterations that are required to 
obtain P(A) = 1. That occurs when the argument of the 
sine function in Eq. (|3.19|1 is tt/2, i.e., 



fc=if^-l 
2 \29 



(3.20) 



For small 9 (or large k) 



k-- 



integer value of — 
46* 



TT 1 

4 7X. 



(3.21) 

Thus after approximately 7r/(4\/A) iterations one has 
certainty of having found the superposition of marked 
states. Classically the number of search operations 
to have this certainty is on the average approximately 
1/(2A) = N/{2M) for N much larger than M. By the 
same reasoning we find P = with twice as many quan- 
tum iterations. Thus by continuing to iterate indefinitely 
we can end with any probability of success. However, if 
we iterate close to the number that gives 100% proba- 
bility of success we have a good approximation to a suc- 
cessful search. 



B. Effect of phase a 



For a general value of a; = (1 — cosa)A, Eq. (|3.12p can 
be written as 



Uk = A cos{k(j) + 9) , 



(3.22) 



where 



2-x 



A=\l'^ — 9 ^ SiictaiiJ -r^. (3.23) 



2-x 



Since P(A) = 1 — u| = 1 — cos {k(j) + 9) , the minimum 
of P(A) is 



2 — (1 — cos a) A 



(3.24) 



In Fig. [S] it is seen that Pmin(A) has at most a linear rise 
as A increases from zero to one. 



IV. SUMMARY 

We have proposed a multi-phase matching for the 
Grover search algorithm, which is an extension of the new 
phase matching proposed in Rcf. [10]. The multi-phase 
matching is characterized by multiple Grover operations 
with two kinds of multi-phases aj and (3j (j = 1, . . . , fc). 
We showed that if we match aj and (3j in accordance 
with the rule aj = —Pk-j+i for a given k we can obtain 
an optimal solution for aj , Pj that gives a success prob- 
ability curve such that it is almost constant and unity 
in a wide range of the fraction of marked states. As an 
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FIG. 8: (Color online) Plots of Pmin(A) for a 
0, 7r/4, Tv/2, 37r/4, and tt - 0.1. 



example we presented an optimal solution obtained for 
k = 6. The solution yields the desired success probabil- 
ity P = 1 to within 0.2% for the fraction of the marked 
states greater than 0.1. This is significant in the sense 
that when the fraction of marked states is greater than 
0.1, we can always with a high degree of confidence find a 
uniform superposition of the marked states by repeating 
the Grover operation just six times. 

To clarify the mechanism of the multi-phase matching 
we studied in detail the one- and two-iteration cases. We 
showed that it is possible to obtain P = 1 exactly for a 
particular fraction A by tuning the phases. This can be 
generalized to having k values of A for which P{X) = 1 
when we go to a fc-iteration scheme. 

One can obtain P = 1 for a given very small A by using 
the original Grover algorithm or the phase-matched ver- 
sion of it. In this case usually a specified large number of 
iterations is required. Further study is needed to obtain 
an efficient algorithm for extremely small A. 
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APPENDIX A: EQUIVALENCE OF TWO 
PHASE-MATCHING SCHEMES 



Li and Li [lO| claim to have generalized the Long phase- 
matching algorithm in order to produce a higher success 



probability. In actual fact the phase matching of Li and 
Li and that of Long et al. @ result in the same success 
probability. We show that in the following. 

Instead of Eqs. (|2.ip and (|2.2p . Long et al. work with 
the operators 



M-l 



U = 1 -{I 



y = /- (l-e*'^)|0^")(0'^"| 



(Al) 



(A2) 



These unitary transformations lead to the Grover opera- 
tor (in the notation of this paper) G{9, (j)) , where 



g ^ / 1 - (1 - e^^ )(l-A) -(1 - e**)e'VA(l-A)\ 
\^-(l-e"^)^A(l-A) [1 - (1 - e'''^)A]e"' J' 

(A3) 

For one operation we calculate the final state 



(A4) 




with 



VT^[1 - (1 - e*^)(l - A) - (1 - e''^)e'^A]. (A5) 



Setting u 
tion 



we obtain (in addition to A = 1) the solu- 



A 



1 cos + 1 

2 sin^d 



(A6) 



Note that the signs of (j) and 9 are the same, unlike the 
opposite signs of the matched phases of Li and Li, i.e., 
/3 = —a. In order that < A < 1 with this phase match- 
ing, 9 varies from 7r/3 to tt. For P(A) = 1 — we 
obtain the expression of Eq. (|2.12[) with a replaced by 
9. Thus the impressive result by Li and Li of a single 
phase-matched Grover operation can also be obtained 
with the earlier-proposed operation of Long et al. How- 
ever, the formulation of Li and Li results in Imu ~ 
when P = —a, whereas Imu ^ when <j) = 9 for the 
operator of Long et al. It should be noted however that 
the remarkable single-operation result was first reported 
by Li and Li [T^]. Although the probabilities are the 
same the amplitudes are not, and Li and Li's formulation 
gives a more straightforward derivation of the probabil- 
ities. (See Sec. IIA.) One can relate the two formula- 
tions by suggesting that instead of the operator acting 
on (Vl — A, >/A)"^ initially, in the case of Long et al. it 
operates on this state multiplied by a phase factor. 
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